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We evaluate the entanglement entropy of exactly solvable Hamiltonians corresponding to general 
families of three-dimensional topological models. We show that the modification to the entropic 
area law due to three-dimensional topological properties is richer than the two-dimensional case. In 
addition to the reduction of the entropy caused by non-zero vacuum expectation value of contractible 
loop operators a new topological invariant appears that increases the entropy if the model consists 
of non-trivially braiding anyons. As a result the three-dimensional topological entanglement entropy 
provides only partial information about the two entropic topological invariants. 
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Introduction: The topological features of the entan¬ 
glement entropy of two-dimensional systems is well un¬ 
derstood m- Tracing a large, smooth region A of a 
gapped system prepared at its ground state gives the en¬ 
tanglement entropy 

SA = a\dA\-bojo, ( 1 ) 

where a and 70 are non-negative constants, | 9 A| is the 
size of the boundary of A and the 0 -th Betti number, Bq, 
counts the number of disjoint boundary components 
A non-zero 70 signals that the model is topological in 
the sense that it supports large loop operators with non¬ 
vanishing vacuum expectation value [ 3 ]. As a conse¬ 
quence it exhibits topological degeneracy and it supports 
anyonic excitations [B]. Interestingly, 70 can be isolated 
by linear combinations of the entanglement entropy cor¬ 
responding to suitable partitions of the system so that all 
the boundary terms from the entropy cancel out mm- 
This entropic feature of the ground state, known as the 
topological entanglement entropy, is one of the main tools 
used to probe topological order in two dimensions. 

The entropic properties of three-dimensional topolog¬ 
ical models have recently attracted considerable inter¬ 
est [THini- It has been shown that, even if three dimen¬ 
sions can support string-like topological excitations, their 
contribution to the topological entanglement entropy is 
equivalent to the one obtained by point-like excitation . 
Furthermore, it has been demonstrated that in three di¬ 
mensions it is possible to have non-universal contribu¬ 
tions to the entropy due to a finite correlation length that 
are constant with the size of the partition [ 5 ]. Neverthe¬ 
less, the constant contribution to the entropy is always 
linear in bo and 61, the 0 -th and 1 -st Betti numbers [8]. 
The latter counts how many topologically inequivalent 
non-contractible loops can exist on the boundary dA, e.g. 
it is zero for the sphere and two for the torus [S] . Intrigu- 
ingly, it has been observed that excitations of some three- 
dimensional models are confining, significantly affecting 
the entropic behaviour of the model [ 9 ]. Nevertheless, 
explicit expressions of the entropy have been found for 
the limited examples of the toric code and the semion 


model [ 5 ] as well as the discrete gauge theories [ 5 ]. 

Here we evaluate the entanglement entropy of Walker- 
Wang models that exhibit rich topological behaviours 
m- These three-dimensional lattice models are analyti¬ 
cally tractable, they have zero correlation length and they 
are fixed points of general families of three-dimensional 
topological models. By evaluating explicitly the geo¬ 
metric entanglement entropy we find the general form 
Sa = q;| 9 A| — 6o7o + ^7i, where 71 is a non-negative 
number. As in two-dimensions the constant 70 deter¬ 
mines the total quantum dimensions of the underlying 
anyonic model of the Walker-Wang model. It also in¬ 
dicates if anyonic excitations can emerge either in the 
bulk of the model or at its boundary [ 9 ]. If 70 > 0 then 
the quantity 71 signals if the anyons of the model braid 
trivially or not. Both 70 and 71 are entropic topological 
invariants that provide information about the topological 
properties of the system. Surprisingly, we find that topo¬ 
logical entanglement entropy does not fully identify the 
topological properties of the model as it fails to distin¬ 
guish between the different types of entropic topological 
invariants. Our work provides qualitative and quantita¬ 
tive investigations of the universal topological character¬ 
istics of three-dimensional fixed point models. 

The model:- Consider an anyonic model with charges 
{!,...,n}, where 1 denotes the vacuum. These charges 
satisfy the fusion rules a x b = where the inte¬ 

gers denote the multiplicity of the fusion channels. 
To each anyon a we can assign a real number, the quan¬ 
tum dimension da, that satisfies dadb = ^ab^c- More¬ 
over, each anyon a has a specific spin, giving rise to a 
complex phase factor 9a when it is rotated counterclock¬ 
wise around itself by 27 r. The fusion and braiding prop¬ 
erties of anyons are described by the F- and i?-matrices, 
respectively, with elements given pictorially by 



Two successive braidings give ^Iso known 

as the monodromy. With the monodromy matrix we can 
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introduce the modularity condition. For modular models 
each anyon a ^ 1 has a monodromy operator 
which is not equal to the identity for at least one charge 
b. In other words a modular model has anyons that braid 
non-trivially with each other. Otherwise the model is 
called non-modular. If all pairs of anyons in a model 
braid trivially with each other then the model is called 
degenerate non-modular (or symmetric). 

We next define 5^^ [HI [T3] , in terms of the following 
quantum mechanical amplitude 


cc _ 
‘^ab — 



1 

V 


E nrj pab,cj^j_ 
0aOb 


a / dadbdj 


( 3 ) 


where V = VEa is the total quantum dimension of 
the model. The amplitudes S^b can be considered as the 
elements of the S-tensor. This tensor is a generalisation 
of the 5-matrix that has elements Sab = Saf^ [H]. The 
modularity condition is ensured if the 5-matrix is uni¬ 
tary nn US]. This allows us to define the modularity 
condition 


^J^'^dacba =Sbl. (4) 

a 

In other words the non-trivial braiding of anyon b with 
some other anyon a causes it to be projected when a 
symmetrisation is applied. 

We now introduce the Walker-Wang models HH. They 
are three-dimensional generalisations of the string-net 
models m- We consider a three-dimensional trivalent 
lattice and adopt a certain anyonic model with n differ¬ 
ent types of charges. We assign an n-dimensional Hilbert 
space at each link, each state corresponding to a Wilson 
line of a certain anyon charge. It is possible to assign a 
Hamiltonian to this lattice that gives rise to a specific, 
possibly degenerate, ground state |$), separated from ex¬ 
cited states by a non-zero energy gap. The Hamiltonian 
energetically penalises the states of the links surrounding 
a certain vertex that do not satisfy the fusion rules. It 
also penalises any non-trivial flux that goes through the 
plaquettes of the lattice m- The ground state of this 
Hamiltonian can be considered as a superposition of ar¬ 
bitrary three-dimensional string-net configurations that 
satisfy the anyonic fusion rules at its vertices. 

Consider the ground state 1$) = El '"^here L 

is a certain string-net configuration appearing with am¬ 
plitude $(L). We can determine these amplitudes in 
the following way m- Apart from the conditions im¬ 
posed in the two-dimensional string-net models that im¬ 
plement topological invariance of string configurations, 
scale invariance and change of basis by employing the F- 
matrices m the three-dimensional Walker-Wang ground 
states satisfy the additional condition 


The configuration (T)_) is a possible string-net in three di¬ 
mensions, but not in two dimensions due to its crossings. 
It has been shown that the behaviour of the point-like 
excitations of the Walker-Wang models strongly depends 
on the modularity condition |^. The modular models 
have all point-like excitations in the bulk confined as 
their non-trivial braiding causes a non-zero string ten¬ 
sion to emerge. On the other hand, the non-modular 
models have deconfined bulk excitations. Nevertheless, 
both models have non-zero entropic topological invari¬ 
ants as we shall see in the following. 

Geometric entanglement entropy: To determine 
the geometric entanglement entropy we evaluate the von- 
Neumann entropy Sa = —Ih/Oa) of the reduced 
density matrix pA of a geometric region A. In two dimen¬ 
sions the entropy Sa is given by Q, where 70 = Inl?^. 
Below we show that in three dimensions the entangle¬ 
ment entropy depends on the topology of the boundary 
dA and on the modularity property of the anyonic model. 

Consider first tracing a simply connected region A out 
of the three-dimensional lattice, which is topologically 
equivalent to a sphere. We take the boundary dA to cross 
the links of the lattice. For convenience we double the 
Hilbert space of the boundary links so they are present in 
both A and its compliment B. We can write the Schmidt 
decomposition of the ground state 1 $) of the system in 
terms of the states of the corresponding regions 

( 6 ) 


where i parametrises the states of all the links of the 


(j)» 


$ 


\ 


A/B/ 


boundary dA, l«il = 1 and y^A/B 
6iA'- The entropy of the reduced density matrix pA = 
trB( I$)($!) corresponds to all possible anyonic configura¬ 
tions of the boundary links. As the system is prepared in 
its ground state these links are subject to the constraint 
that they all fuse to the vacuum [16]. The probability of 
having a certain link configuration i of anyonic states at 
the boundary dA is given by P{i) — |ajp. To evaluate 
P{i) we first note that the probability of measuring a 
given charge a at a certain link is given by Pa = da/T>^. 
To impose the fusion constraint at the boundary we in¬ 
troduce the conditional probability of fusing N anyons 
to the vacuum is given by P{ai x 02 x ... x —>■ 1 ) = 

^a/iUi dat)- Here is the to¬ 

tal multiplicity associated to the fusion of the N anyons. 
Therefore the probability that the configuration i occurs 
at the boundary with vacuum total charge is given by 
P{i) = Pa^...Paj.,P{ai X ... X On ^ 1)/Hi [17]. In the 
case where there are several disjoint boundaries compris¬ 
ing dA then each boundary component carries vacuum 
total charge. Hence, the total probability P{i) is the 
product of the individual probabilities giving finally 



( 5 ) 


2 _ -Et riisi dai 
X)2iN-bo) 


( 7 ) 
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The eigenvalues of the reduced density matrix pA are 
given by = laip/A/"^, where cr = with /x 

parametrising the multiplicities of the fusion channels. 
Hence, the entanglement entropy is given by 

= ( 8 ) 

k 


This is the same behaviour as in the two-dimensional 
case 0. where \dA\ = TV is the number of links on the 
boundary of H [3]. 

Consider now the case where the boundary dA is topo¬ 
logically equivalent to a torus, as shown in Fig. (Left). 
Compared to the entropic behaviour of simply connected 
regions the torus can get additional contributions emerg¬ 
ing when loops that have support in A and loops that 
have support in B are braided 0 . To facilitate the calcu¬ 
lation of the entropy Sa we bring the state of the system 
to a suitable form by applications of loop operations that 
have support exclusively in A or in B. These operations 
do not change the value of Sa- We can thus bring the 
braided loops that belong to different regions in the form 
shown in Fig. (Left). There, a/b is the total anyonic 
charge of the non-contractible loops that have support in 
A/B, respectively, and c is the exchanged anyon, as dic¬ 
tated by the 5-tensor. The charges of a, b and c cannot 
be changed by non-boundary operations. Thus, accord¬ 
ing to ([^ the ground state of the system is entropically 
equivalent to 


1 $) 


— ^ ai^S2b 




( 9 ) 


where if, denotes the boundary configuration with a 


charge c crossing the boundary and 
are basis states satisfying ( 


A,a,c 


and 




AjB^ajb-jC 


$ 


Aj B.a' jb' ,c' 


jb'^c,c' ■ Notc that |^) is properly normalised 
as 


= (10) 


due to the unitarity of the F-matrix [TS]. We next 
show that the entanglement entropy of a degenerate non- 
modular model corresponding to a toroidal boundary is 
the same as for the simply connected region before we 
turn to the modular case. 

To calculate the entanglement entropy of the degener¬ 
ate non-modular case we first determine the 5-tensor. As 
these models have trivial monodromy between all anyons 
S'/b — bcidadb/T). Hence, the state ([^ takes the form 


1^) = ^ 


ii,a,6 


^ A.a.l 


- 6 , 6,1 


( 11 ) 


By introducing the orthonormal states 




y- da 

v 


^x4,a,l 


and 




— ik 





FIG. 1: (Left) A toroidal region traced out of the lattice. 
A possible string-net configuration is depicted with anyon a 
threading through the inside of the torus and b threads around 
it. Anyon c is a possible exchanged anyon between the two 
loops according to the 5-tensor. (Right) Example of a possi¬ 
ble configuration of regions A, B, C and D that extract the 
topological entanglement entropy in three dimensions. 


rewrite the ground state 0 in the same form as §. 
Hence, for a degenerate non-modular model the entropy 
of a toroidal region is the same as the entropy of a simply 
connected region given in 

We now consider the case of modular models. To write 
the ground state 0 of the modular case in a more con¬ 
venient form we introduce the following states 

1 CO (f,*c 

Z-^a ^ab ^A,ac 

to the relation 


These states are orthonormal due 


^(5=J*5“,, = dJbb', 


( 12 ) 


which can be shown by employing the modularity con¬ 
dition 0 [nmi]. Then the ground state of the system 
becomes 


ic^b,c 


$ 


A,be 


$ 


B,bc 


(13) 


We can proceed in the evaluation of the entropy in a 
similar fashion as we did for the simply connected case. 
Note though that the condition of having the total any¬ 
onic charge at the surface dA being the vacuum, 1, is 
now replaced to c due to the braiding of the a and b 
anyons. The associated probability P{ic) = for 

the boundary configuration can be written as P{ic) = 
Pai---Pai-rP{ai X ... X Oat —>■ c)/Pc- Combined with the 
rest of the normalisation factors of state (13) we obtain 


the eigenvalues of the reduced density matrix pA to be 


A. = 




A/f 




(14) 


We can directly extend these probabilities to the case 
of higher genus regions with boundary Betti numbers bg 
and bi by introducing the 5-tensor in 0 multiple times, 
giving finally the entropy 


Sa = -N 


^ di 4 


^ X>2 

k 




(15) 


we can 
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To derive this result we have used the relation 
= 1 , which follows from ( 10 ) and ( 12 ). For 


two-dimensional surfaces, such as the boundary ijA, the 
Betti numbers are related by x = 2 &o ~ ) where x is the 

Euler characteristic. Hence, the topological dependence 
of the entropy for modular models is proportional to the 
Euler characteristic. 

We now consider a family of models with entropic be¬ 
haviour that ranges between the modular and the degen¬ 
erate non-modular cases. Take two anyon models, a de¬ 
generate non-modular one ^Idnm and a modular one 
with quantum dimensions I?dnm and I’m, respectively. 
We can construct a new anyon model as the direct prod¬ 
uct A = -Adnm X The charges of such a model are 
given by a = {oi, 02}, oi S ^dnm, 02 S ^m, with = 

total quantum 

dimension V = PdnmI^m- As such = ‘Salbi‘^Zb2- 
the two anyon models are non-interacting the entropy of 
A can be written as the sum of the entropies of each 
component, 5'2i(A) = S'a(Adnm) + Then the 

entanglement entropy is given by 


Sa = -N 




X>2 

kGA 




M- 


( 16 ) 


In this way we can construct anyonic models with arbi¬ 
trary entropic behaviours, ranging between the modular 
(2 ?dnm = 1) and the degenerate non-modular (I?m = 1) 
cases. 

Entropic topological invariants:- The entangle¬ 


ment entropy ( 16 ) provides the general form of the en¬ 
tropic topological invariants, 70 = Inl?^ and 71 = Inl?^. 
The first invariant, 70, uniquely determines if the model 
is topological, 70 0 , or not, 70 = 0 . Similarly to the 

two-dimensional case it indicates if loop operators, cor¬ 
responding to arbitrarily large but contractible loops can 
have a non-zero expectation value. In accordance to the 
anyonic fusion rules it reduces the entanglement entropy 
generated by the superposition of the string-net states. 
The second invariant, 71, is always less than or equal to 
7o. It quantifies the total quantum dimensions of anyons 
in the model which braid non-trivially. If 71 = 0 then 
the model is degenerate non-modular, and all its anyons 
braid trivially with each other. These anyons can exist 
in the bulk of the model as deconfined excitations. If 
7o > 7i the model is non-modular supporting both con¬ 
fined and deconfined excitations [ 3 ]. If 70 = 71 >0 then 
all excitations of the model are confined and the model 
is modular. Note that excitations that are confined in 
the bulk can appear deconfined at the physical bound¬ 
ary of the system nms]. Topological degeneracy can be 
manifested by imposing periodic conditions to the three- 
dimensional bulk or its two-dimensional boundary [ 9 ]. 

The natural question is what information is gained 
from the three-dimensional generalisation of the topolog¬ 
ical entanglement entropy. Similarly to two-dimensional 


models it is possible to partition the Walker-Wang lattice 
in four regions. A, B, C and D, such that a combination 
of their entanglement entropies gives rise to topological 
invariant quantity. An example of such partition is shown 
in Fig. (Right). To define a topological invariant of en¬ 
tropies we first introduce the quantity [ 5 ] 


S[X(A, B, C)] =Xa + XbFXc-Xab - Xac - Xbc 
+ Xabc, ( 17 ) 


where X is some property of the model. For X = S 
all surface contributions of (161 cancel out apart from 
a possible intersection of the four regions [8]. Hence, 
X)[S{A, B,C)] is a topological invariant for a choice of 
regions that do not have such intersections. In this 
case D[S{A, B,C)] is called the topological entangle¬ 
ment entropy 7. If we now calculate the same quan¬ 
tity for the Euler characteristic, X; of the boundary 
we find D[x( 9 (A, B, C))] = x(A Ci B H C H D) for any 
choice of four regions m- Demanding topological in¬ 
variance of D[S{A, B,C)] gives 'D[x{d{A,B,C))\ = 0 . 
From ( [1^ we obtain the topological entanglement en¬ 
tropy 7 = I)[ 6 o( 9 (A, H, C))] ln(DQj,j]y[ that can identify 
the presence of deconfined excitations in the bulk. This 
analysis also shows that in three dimensions it is not 
always possible to isolate the topological invariant quan¬ 
tities 7o and 71 by combining the entanglement entropies 
of a fixed partition. 


Conclusions:- We have seen that the entanglement 
entropy of the three-dimensional Walker-Wang models 
gives two distinct entropic topological invariants, 70 and 
7i. While 7o identifies if the model is topological, 71 
reveals information about the braiding properties of its 
underlying anyonic model. The latter corresponds to an 
increase in the entropy as tracing the region around a 
torus erases information about the anyonic charge of the 
threaded loops that could affect non-trivially the state 
of the system inside the torus. Hence, the geometric 
entanglement entropy in three dimensions provides more 
information about the topological order of the system 
than its two-dimensional counterpart. We can isolate 
7o by comparing the entropy of two different partitions, 
with different number of disjoint boundary components, 
bo, but the same bi and | 9 A|. Similarly we can determine 
7i. Evaluating the entanglement entropy for the general 
non-modular case is a complex problem due to the lack 
of structure of these models. We leave this problem to 
future investigation. 
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Supplementary Material 


A. Probability of boundary configurations:-There 
are two equivalent ways to assign a probability to a given 
boundary configuration on a region A of a string-net con¬ 
figuration in Walker-Wang models. The first is to use 
unitary operations with support on either A or the com¬ 
pliment B. One can then form a canonical configura¬ 
tion with string-nets restricted to the plaquettes crossing 
the boundary dA. Then the diagram calculus can be 
utilised to calculate the amplitude of each configuration 
and hence the probability [ 3 ]. The second approach is 
to assign a conditional probability to the anyon config¬ 
uration making use of the quantum dimensions of the 
model [ 2 ]. This is the methodology we utilise here as it 
is computationally easier . 

In order to calculate the probability of finding an anyon 
charge we make use of the iS-matrix. The probabil¬ 
ity of finding an anyon of charge a in a particular re¬ 
gion is given by Pa = = da/D"^ [ 2 ]. This result 

can be understood as the probability amplitude for the 
vacuum braiding trivially with the charge a. Due to 
charge conservation across the boundary the probabil¬ 
ity of each anyon on the boundary is not independent 
but we require that all anyons fuse to a total charge 
c. For a boundary consisting of N charges this con¬ 
straint is applied by defining the conditional probability 
P(ai X ... X On c). This probability can be evaluated 
diagrammatically making use of isotopy invariance |20j 

and the identity = \fS; ffH] 


P(aia2 —f ji) — P{aia2 —>■ ji )^—^— 

^ai^a 2 




aia2 Ji 

^ai ^02 


(Al) 


We now calculate P(ai x ... x oat —>■ c) as below 


P(ai X ... X Oat —1 c) 

= P(aia2 —f ii)P(a3ii j2)---P{o-NjN-2 c) 
_ J^gdc 

■prA^ 7 ’ 

llj^l ^ 0.1 


(A 2 ) 


where Mf, = can verify the 

consistency equation 

Pa--PbP{ai X ■■■ xpN ^ c) = Pc. (A 3 ) 


The normalised probability for the boundary configura¬ 
tion with a given charge c across the boundary is given 
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by 


P(iJ = X ... X ^ e) 

Pc 

B. S-tensor Properties:- The 5-matrix has been 
extensively studied throughout the literature [Til [T5] . 
Here we outline the generalisation to the 5-tensor and 
we demonstrate some of its properties. The 5-tensor can 
be defined diagrammatically or in terms of the parame¬ 
ters of the model as follows m 



1 

V 


E jsrj pab,cJli_ 
“''d OaOb 


a/ dadbdj 


(Bl) 


One useful property of the 5-tensor is the following state¬ 
ment 




(B2) 


This follows from Eqn. (Bl) noting the unitarity of the 
F-matrices J2c(^ab'j)*^abj' — ffS]- In Ibe remain¬ 
der of the section we make use of the diagram calculus 
in order to derive properties of the 5-tensor. 

We begin our discussion by defining the fusion Hilbert 
space and the relevant normalisation scheme [131 [El ED] ■ 
For the fusion process a x 6 —>■ c we can define a fusion 
Hilbert space with dimension given by the fusion 
multiplicities =dimiy“*’. We define an orthonormal 
basis for such a space with states |a, &;c, sat¬ 

isfying (o', 6 '; c',/x'l a,b]c,p.) = 5aa'5bb’bcc'5^^'. We can 
represent these states and their duals in terms of the di¬ 
agram calculus 


a b c 

\a,b;c,p)= , (a, 6 ;c,/r|= . (B3) 

c a b 

The inner product is dehned diagrammatically by stack¬ 
ing vertices vertically and connecting open edges. To 
preserve isotopy invariance and normalise the inner prod¬ 
ucts we introduce a constant depending on the quantum 
dimensions of the edge charges 


We now represent the 5-tensor as an operator 5c on the 
fusion space as defined in [13] 


Sc\b,c-b,fj.) 

X 



(B 6 ) 


Utilising the diagram calculus and normalisation conven¬ 
tions dehned above we can form the matrix elements of 
the 5-tensor operator 

Sab = \/~dc{a,c-,a,p,\Sc\b,c-,b,p) . (B7) 

The above discussion of the 5-tensor holds true for any 
anyon model. In the following we restrict our discussion 
to modular anyon models to show that for such models 

Y.{si,,rsib = dc5b,b'. (B 8 ) 

a 

As discussed in the main text an anyon model is described 
as being modular when for each non-vacuum charge a in 
the model there exists a charge b such that the mon- 
odromy 7 ^ lab- This statement can be equiva¬ 

lently formulated in the diagram calculus as the modular 
trap identity [13 E] 


1 ^ yx 

a I 



X X 


w 


X X 


(B9) 


Using this identity it is possible to show that the operator 
5c is unitary for modular models 


5l5c|&,c;5,^) qJJ/ 


x,V 


/—b 


b c 

= ^ =\b,c;b,n). 

b 


(BIO) 


We then prove (B 8 ) by invoking the unitary property of 
5c and (jBTl) 


a , 6 ; c , fj, \ a,b; c, fi) = 


^aa' bbb'bcc'd^ 

y/dadbdc 



— daa'dbb'dcc'd^fj,' ■ 

We also introduce the identity operator 


(B4) 


a b 


1= {a,b-c,fi\= % ■ 

a.b,c,M a,b,c,ii * “ 


a b 

(B5) 


dcdbh' 

=dc (&',c;6',/i| b,c;b,p.) 

=dc (&', c; 6', fi\SlSc \b, c; 6, fi) 

= Y dc{b',c;b',fj,\Sl\a,c';a,n'){a,c';a,p,'\Sc\b,c;b,p,) 

a,c' 

= Ydc {b', c; b', p.\ si \a, c; a, p) {a, c; a, p\ 5c \b, c; b, p) 

a 

^Yisib'Y^:,. (Bii) 

a 

C. Isolating topological components of the en¬ 
tanglement entropy:- As for two-dimensional models 
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[3 la, it is possible to partition the Walker-Wang lat¬ 
tice into four regions A, B,C and D from which a com¬ 
bination of entropies can be evaluated to formulate a 
topologically invariant quantity [71 [5]. To define a topo¬ 
logically invariant quantity we introduce the quantity 
T)[X{A^B,Cy\ as in [5] 

D[XiA, B, C)] =Xa+Xb+Xc-Xab - Xac - Xbc 
+ Xabc, (Cl) 


We utilise the above construction to define the boundary 
of the regions A, B, C and their composites with respect 
to union and intersections of the other regions 

d{ABC) = {ABC) n D, 
d{AB) = {AB) n {CD), 

dA = Af^{BCD). (C2) 


where X is a property of the regions A, B, C and D. 
Work by Grover et al. [8] conjectured that for X = S' all 
surface contributions from the geometric entanglement 
entropies cancel when the quantity S)[x(9(A, B, C))] = 0, 
where x is the Euler characteristic. In the following we 
show that under minimal constraints, surface contribu¬ 
tions and terms proportional to the Euler characteristic 
vanish identically in equation (Cl). Nevertheless, our 
workings show that the vanishing of the Euler character¬ 
istic does not imply the vanishing of surface contributions 
as claimed in [H]. 

Consider the lattice in the continuum limit as a closed 
subset of IR^ partitioned into four finite regions A, B, C 
and D with non-intersecting volumes where A, B, C are 
embedded within D. We define the condition of non¬ 
intersecting volumes by the statement, two volumes I 
and J are non-intersecting if and only if / n J = dlC\ dJ 
where dl is the boundary of I. By convention we adopt 
the following notation for composite regions IJ = IU J. 


The boundaries of the other regions can be formulated 
likewise. 

Measures of the surface area and Euler characteristic 
both obey the so called inclusion-exclusion principle. A 
property of a set X obeys the inclusion-exclusion princi¬ 
ple when X{AB) = X{A) + X{B) — X{AD B). Making 
use of this property one can easily verify 

T)[x{d{A,B,C))]= X{AnBnCnD). (C3) 

When X is taken to be the surface area or Euler charac¬ 
teristic, Area(0) = x(0) = 0, hence for AflBnCdD = 0, 
'D[X{d{A, B,C))] = 0. This result demonstrates that 
such a combination of boundary areas or Euler character¬ 
istic will always vanish for the above construction when 
there is no simultaneous intersection of all four regions. 
Furthermore such a quantity is necessarily invariant un¬ 
der smooth deformations of the regions complimentary 
to the discussion in [8]. 



